Creep motion of a model frictional system 
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We report on the dynamics of a model frictional system submitted to minute external perturba- 
tions. The system consists of a chain of sliders connected through elastic springs that rest on an 
incline. By introducing cyclic expansions and contractions of the springs we observe a reptation of 
the chain. We account for the average reptation velocity theoretically. The velocity of small systems 
exhibits a series of plateaus as a function of the incline angle. Due to elastic effects, there exists 
a critical amplitude below which the reptation is expected to cease. However, rather than a full 
stop of the creep, we observe in numerical simulations a transition between a continuous-creep and 
an irregular-creep regime when the critical amplitude is approached. The latter transition is rem- 
iniscent of the transition between the continuous and the irregular compaction of granular matter 
submitted to periodic temperature changes. 

PACS numbers: 45.05.+X General theory of classical mechanics of discrete systems, 45.70.-n Granular sys- 
tems, 46.55.-|-d Tribology and mechanical contacts, 65. 40. Do Thermal expansion; thermomechanical effects. 



I. INTRODUCTION. 

Granular materials are collections of macroscopic par- 
ticles (grains) that interact via dissipative forces. As a 
result, external excitation is necessary to promote the 
motion of the grains. In general, it is assumed that 
in absence of mechanical perturbations, such materials 
will eventually reach a mechanically static configuration 
and remain at rest. In particular, thermal agitation will 
not suffice to induce rearrangements and, for this rea- 
son, granular matter is said to be athermal [IHSl ■ This is 
however a simplified view that assumes an idealized situ- 
ation in which all perturbation (temperature variations, 
humidity changes, mechanical noise, etc.), even minute, 
can be suppressed. 

In practice, small temperature variations can actually 
trigger small but measurable rearrangements of the struc- 
ture. Uncontrolled thermal dilations have been reported 
to generate stress fluctuations large enough to hinder re- 
producible measurements of the stress field inside a gran- 
ular pile [J, |5| . They were even suspected to be the driv- 
ing factor leading to large-scale 'static avalanches' @. 
When accumulating, such minute perturbations can in- 
duce an irreversible evolution of the system. Several re- 
cent studies indeed showed that temperature cycles, even 
of small amplitude, can induce the slow compaction (a 
creep) of dry granular materials (7l-[l3|. Moreover, a tran- 
sition between a continuous-flow and an intermittent-flow 
regime, observed when the amplitude of the temperature 
cycles is decreased [l0l - [l3| . remains unexplained. Even 
if the transition is thought to be due to finite size effects 
{i.e. the finite number of grains in the diameter of the 
tube), the mechanisms, in particular the role played by 
the confining walls, are still under debate. 

In the same manner, minute temperature changes have 
been identified, more than a century ago, to induce the 
creep of solids in frictional contacts. H. Moseley reported 



first the descent, driven by temperature variations, of 
lead plates covering the south side of the choir of Bris- 
tol cathedral and mentionncd that the same mechanism 
could be responsible for the motion of glaciers [l^l- A 
macroscopic model, considering only the relative dila- 
tions of the solids in contact, predicts that the creep ve- 
locity is prop ortional to the amplitude of the temperature 
changes jl4l ITsI and, thus, that any temperature varia- 
tion, whatever its amplitude, leads to a motion. However, 
as noticed by H. Bouasse, such a description does not take 
into account the elastic effects which necessarily play a 
role for large systems • When elastic effects come into 
play, the temperature variations are sufficient to induce 
the motion only if their amplitude is larger than a criti- 
cal value which, in particular, increases with the size of 
the system [l6j . Accordingly, for a given amplitude of 
the temperature variations, a solid is observed to descent 
along an incline only if the tilt angle is large enough [17 1 . 



In the present article, we propose the study of a model 
mimicking the reptation along an incline of a solid sub- 
jected to temperature variations. We remark that the 
aforementioned models disregard the fact that the sur- 
faces, even if nominally flat, are in real contact in a finite 
number of localized regions. However, due to roughness, 
the real contact between the surfaces reduces to a large, 
but flnite, number of microcontacts, themselves belong- 
ing to a small number of mesoscopic, coherent, contact 
regions [isj . Such coherent regions are composed of a 
large number of microcontacts, so that their contact with 
the substrate obeys the Amonton-Coulomb law for static 
friction. We thus consider a set of sliders, connected 
by springs, in frictional contact with an incline and sub- 
jected to temperature changes. We report and discuss an 
extensive study of the associated dynamics. We point out 
that, interestingly, the conclusions help in understanding 
the observations reported for a granular material in a 
tube. 
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II. THE MODEL 
A. Description 

The system (total mass M) consists of a series of N 
identical sliders (zero length, mass m = M/N) connected 
by {N — 1) identical springs (massless, stiffness k). The 
sliders are in frictional contact with an incline which 
makes an angle, a, with the horizontal (Fig. [T]). They 
are subjected to the forces due to the springs, to their 
own weight mg (where g denotes the acceleration due to 
gravity) and to the reaction force from the incline (in- 
cluding the frictional force). 




FIG. 1. Sketch of the model system. 



The aim of the study is to account for the dynamics 
of the system induced by temperature changes. To do 
so, we consider that the natural length of the springs, /, 
depends linearly on the temperature T according to: 

1{T) ^ + n{T - To)] (1) 

where k stands for the thermal expansion coefficient of 
the slider material and Iq for the natural length of the 
springs at Tq. In accordance, the total length of the sys- 
tem at To, in absence of internal stress, is i = {N — 1) Iq. 
We further assume, for the sake of simplicity, that the 
stiffness k does not depend on the temperature and that 
the substrate does not dilate. 

^From now on, the x-axis is oriented downwards, a 
is positive and the first slider is the upper one (Fig. [1]). 
Thus, at the temperature T, the force due to the spring, 
exerted by the slider {n + 1) on the slider n is: 

fn+l^n = -k [Xn - X^+l + l(T)] (2) 

where a;„ denotes the position of the slider n on the in- 
cline. 

At rest, the contact between the sliders and the incline 
is characterized by the static frictional coefficient fis,n 
such that the slider n, initially at rest, starts moving if: 

l/ri+i^n + fn~i^n + mg sina| > ^jis,n mg cosa. (3) 

We remark that, at the boundaries, /o->.i = (n = 1) 
and Jn+i^n = (n = N). Note that, due to the spa- 
tial heterogeneity of the incline surface, the static fric- 
tional coefficient might take, at random, different values, 
fis,n, for the different sliders which lie at different posi- 
tions, Xn, on the incline [l^. In accordance, a slider that 
moves and stops at a different position might be then 
associated with a different value of fJ.s,n- By contrast, we 



assume that, when the slider is in motion, the frictional 
contact is characterized by a single, constant, dynamical 
frictional coefficient, /id. Indeed, while in motion, the 
slider explores the incline and is thus less sensitive to the 
details of the surface properties. The constant /id quan- 
tifies the average rate of energy dissipation and the slider 
n in motion is subjected to the associated dynamical fric- 
tional force: 

fd,n = -fJ-dmgcosaS{xn) (4) 

where i„ denotes the velocity and S the sign function 
[S{u) = 1 if u > and S{u) = -1 if u < 0]. Wc further 
assume, in agreement with standard observations, that 
fid < fJ's,n (Vn), i.e., the dynamical frictional coefficient 
is smaller than the static one. 

The mechanical system is submitted to temperature 
changes. Some sliders, initially at rest, start moving 
when the length 1{T) exceeds a value such that the con- 
dition ([3|) is satisfied for, at least, one slider. Submitted 
to the forces due to the springs, to their own weight and 
to the reaction force from the incline (including the fric- 
tional force), one or more sliders move and come back 
to rest. Such event constitutes, by definition, a ^'me- 
chanical rearrangement^ of the system. At this point, it 
is interesting to consider two characteristic times. On 
the one hand, due to its thermal inertia, the system ex- 
hibits a thermal characteristic time rth which limits the 
dynamics of the temperature changes and, thus, of the 
thermal dilations. In practice, for a macroscopic system 
whose typical size L ranges from a few millimeters to 
centimeters, Tth ranges from seconds to hours. On the 
other hand, the dynamics of the mechanical system is 
characterized by the typical time scale Tdyn = \/ mjk. 
In practice, one can consider instead that Tdyn L/cs 
where Cg stands of the speed of sound in the material the 
macroscopic solid is made of. For a typical size L ranging 
from a few millimeters to centimeters and usual values 
of Cg (about a few kilometers per second), we estimate 
Tdyn ^ 10^^ — 10~^ s, thus much smaller than rth. As a 
consequence, we consider that /(T) is constant during the 
mechanical rearrangements which are, in practice, much 
faster than the temperature variations. 



B. General system of equations 

1. Mechanical system 

First, wc write the differential equation governing the 
position a;n of the slider n, when the latter is in motion. 
Introducing the thermal dilation 9 = k(T — Tq) and the 
dimensionless variables t = 1 1 Tdyn and x = x/{gT^ ) 
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(with Tdyn = \/m/k), we get: 

xi = - [xi - X2 + {1 + 6) la] 

— /id >S'(ii) cosa + sina (n = 1) (5) 

Xn — [2 {x^j^x Xji—\Y\ 

— /id 5'(a;„) cosa + sina {n^l,N) (6) 

XN = ~[XN ~ Xpf-l - (1 + 0) Iq] 

— /id '5'(5;Af) cosa + sina {n ~ N) (7) 

where, in accordance, we defined Iq = ^o/(ff T'dyn)- 

Second, we remind that the shder n, if at rest, starts 
moving if the condition (jSj is satisfied. Thus, using the 
dimensionless variables, we can write: 

\x2 ~ xi — (1 + 9) Iq + sinaj > /is,i cos a (8) 
\xn+i + Xn-1 - 2 x„ + sin a\ > /is,„ cos a (9) 
|i:Ar_i — xn + {1 + 9) Iq + sin a] > ^s.n cos a (10) 

where the conditions ([8]) and ^T0\\ apply for the slider 1 
and respectively and the condition ^ for any other 
slider n 1,N. 



2. Frictional contact 

We remind that the dynamical frictional contact be- 
tween the sliders and the incline is characterized by a 
single coefficient, //d- By contrast, each time a slider 
in motion comes back to rest, a new value of the static 
frictional coefficient /Xg is drawn from a probability dis- 
tribution p(/is)- In order to evaluate if the fluctuations in 
/is are relevant for explaining the qualitative behaviour 
of the system and to be able to discuss the results analyt- 
ically, we assume that p(/is) is a Gaussian distribution, 
namely 



p(Ms) 



1 



: exp 



(Ms - Hs)^ 

2-^ 



(11) 



where /Jg denotes the average value and ct^ the width 
of the distribution. We point out that, the distribution 
p{Hs) does not, strictly speaking, satisfy the condition 
pifJ-s) = 0, V/^s < fJ-d (i-C, the static friction coefficient 
can be smaller than the dynamic one). However, we re- 
strict our study to the case <^ (/Ig ^ A*d), so that the 
latter condition is reasonably satisfied in practice. 



3. Temperature variations 

The system is driven by temperature variations which 
induce changes in the natural length of the springs. We 
will either consider that the temperature oscillates with 
the period 2Tth between two well-defined values such 
that the dilation oscillates periodically between ~Ag and 
+Ae, or that the temperature T of the system fluctuates 



around the temperature Tq with a typical amplitude AT 
such that 9 fluctuates around zero with a typical ampli- 
tude ag ^ K AT. In this latter case, in order to be able 
to discuss the results analytically, we shall assume that, 
at a series of equally-spaced timesteps tq = qrth (g G N), 
a value 9q of the dilation is drawn from the Gaussian 
distribution: 



: exp 



9^ 
2<jI 



(12) 



where ag accounts for the typical amplitude of the ther- 
mal dilations. We then assume that 9 



9q for t e 



[tq,tq+l[. 



C. Numerical method 

Consider that, when the time t reaches tq+i, the sys- 
tem already experienced k mechanical rearrangements 
such that all the sliders are resting at the positions 
x'^ associated with the static frictional coefficients /ig„ 
(ri S [1, A'^]). At t ~ iq+i; a new value of 6 is drawn 
at random according to the distribution ip{9). The sub- 
sequent evolution of the system is obtained as follows. 

From the conditions ([5|)- ([TU1) . one calculates the first 
critical value 9^'^-^ which leads to the destabilization of 
at least one slider (in general, |6'|^"''^| < l^g+ij). Then, 
for 6 = 9^'^^, the equations ©-([T]) arc integrated nu- 
merically using an integration time step At <^ Tdyn (we 
use a standard velocity Verlet integrator [2l|). In order 
to take into account that the motion of one slider can 
destabilize its neighbours, the procedure checks if any 
condition ([5))- PII)) for the sliders at rest is fulfilled within 
each timcstcp At. If so, the motion of the correspond- 
ing slidcr(s) is also accounted for through (O-©. In the 
same way, the procedure checks at each timestep if any of 
the sliders in motion comes to rest. If so, a new value of 
the associated static frictional coefficient is drawn at ran- 
dom according to the distribution p(/is) [Eq. ([TT|) ]. We 
consider that the mechanical rearrangement ends when 
all the sliders are back to rest. All the sliders are then 
resting at the new positions i^"*"^, associated with a new 
set of static frictional coefficients /ign^- 

The procedure is iterated by finding, in the same way, 
the next critical value of the dilation 9c which leads, ac- 
cording to the conditions ([5))- ((TU)) . to the next mechanical 
rearrangement. For each 0c, a new static state of the sys- 
tem is found. The procedure stops when the next critical 
dilation 0c is beyond 9q+i. The system then experienced 
k' rearrangements. The are then the positions of the 
sliders at t = i<j+2- 

The long-term behavior of the system is assessed by 
iterating the whole procedure, either making oscillate 
between ~Ag and +Ag or drawing at random the next 
value of 9 according to the distribution ip{0). 
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III. RESULTS 

A. The minimal system: 2 sliders 

Consider first the system consisting of two sliders 
(N — 2) connected by one spring. The system is very 
similar to that already studied in Ref. [201 fo'" the case 
of a horizontal substrate but, now, the system lies on an 
incline which makes a finite angle a with the horizontal. 



1. Numerical results 

In order to account for the creep of the system along 
the slope we consider the position of the center of mass 
xg = jf ^n=i (here = 2) as a function of time at 
the timesteps tq [Fig. ^] for Gaussian variations of the 
temperature (Eg. [T^. We observe that, for a sufficiently 
large ag, the center of mass moves downwards with, in 
average, the constant velocity < vq >. 



a wide range of a, < vg > increases only slightly with 
a (plateau) until a second critical value a+ is reached. 
Above , < Vg > drastically increases as a approches 
the critical angle of avalanche, ac, defined by tan(Q;c) = 
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FIG. 2. Creep motion of 2 sliders — Top: The dilation 
6 evolves randomly according to the Gaussian distribution 
(|12|l . Bottom: In response to the dilations, the whole systems 
creeps along the incline. We report the displacements Aii 
and Axg of the sliders, as well as the displacement Axg of 
the center of mass with respect to their initial positions. We 
observe that the lower slider moves in response to dilations 
whereas the upper slider moves in response to contractions. 
Note that both sliders move systematically downwards. Inset: 
Long-time behavior - Axg increases linearly with time t, the 
system is creeping at constant velocity (Zo = 10^, ^d=0.5, 
7I,=0.6, a^, = 0.01, ae =0.1 and tana = 0.25). 

For a given amplitude of the temperature variations 
ue, < Vg > 8& & function of a (Fig. [3]), we observe that 
the average creep velocity of the center of mass depends 
only slightly on the slope of the incline over a wide range. 
Indeed, for small a, < vg >; which equals zero for a = 
0, increases rapidly with a and reaches a plateau for a 
typical value that we denote a~ . Then, above a~, over 



FIG. 3. Velocity < vg > vs. angle a — The average 
velocity of the center of mass, < va > depends only slightly 
on a (plateau) over a wide range of the incline angle typically 
from to at- Below aj, < va > vanishes for vanishing a. 
Above at, < va > drastically increases when a is increased 
and diverges for a — ac {lo ~ 10^, /id=0.5, 11^=0.6, = 0.01 
and ae = 0.01). 



2. Analysis 

The results presented above can be qualitatively un- 
derstood by simple arguments. 

First, for a = 0, by symmetry, we expect < vg>= 0. 
This case has been extensively studied in Ref. [201 ■ The 
system ages and the center of mass experiences a sub- 
diffusive motion. Second, for a > ac — arctan(/id), if the 
system starts sliding as a whole, the x-component of the 
weight exceeds the dynamic friction force and the system 
accelerates. A finite average velocity is not defined in this 
case. 

The most striking result, i.e., the plateau observed over 
a wide range of the incline angle (a~ ^ o; < ct^), can 
be easily understood by considering the mechanical sta- 
bility of each of the sliders. Indeed, assuming that the 
values of the static friction coefficients /Xs,i and ^s,2 are 
close to each other, one can see that, because of the force 
mg sin(a) due to the gravity, the lower (resp. upper) 
slider moves downwards when the system dilates (resp. 
contracts). As a consequence, ifOq+i > 9q (the system di- 
lates by A9 = Oq+i — 6q > 0), the resulting displacement 
of the center of mass is Axg = \^X2 — i^lviAO. In the 
same way, if Oq+i < Oq (the system contracts by A6' < 0), 
Aic = ^Axi ~ — iZoA6'. Thus, whatever the sign of A6' 
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(a contraction or a dilation), Axg ~ i/o|A0|. Note that 
Aic does not depend on the angle a. Thus, when the 
system is subjected to periodic dilations from ±Ag to 
T^e, the center of mass is displaced by Axq ~ ^lo{2Ag) 
during Tth such that the average velocity of the system is 
given by: 



Tth < €'G > - loAe 



(13) 



In the case of Gaussian fluctuations [Eq. ([T^]. taking into 
account the distribution 'v^^ estimate the average 

velocity of the center of mass for a~ ^ a < : 



Tth < wg > ~ —^=loO'e- 



(14) 



The latter theoretical predictions ([T3| and (|T4l) are in 
good agreement with the numerical results for large val- 
ues of Ag and ag in a finite range of a (Fig. [3]). 

Let us focus on the critical angles a~ and which 
limit the plateau. On the one hand, the plateau ve- 
locity is reached if the angle a is large enough {a > 
a~) for the force due to the gravity to insure that 
the upper slider remains stable when the system di- 
lates and, conversely, that the lower slider remains sta- 
ble when the system contracts. For instance, for a di- 
lation, the lower slider (the slider 2) destabilizes first if 
/i_>.2 -I- mg sin a > ^s,2 mg cos a [Eq. (fTO|) ] and /2->i + 
mg sina > —fXs^img cosa [Eq. ([5])]. The condition is 
thus, by summing the two inequalities, tana > ^(Ais,2 — 
/is,i). Taking into account the probability distribution 
p(/^s), wc deduce that the latter condition is generally 
fulfilled if tana > tan a^T = ct^./\/2, which defines a^T 
(Fig.©. 

On the other hand, the velocity of the system is larger 
than the plateau velocity if the motion of one slider 
can induce the motion of the second one (a > o.t)- 
The angle a+ can be estimated by considering, for in- 
stance, a dilation such that the slider 2 starts moving 
for fi-f2 + mg sina > Hs,2'mg cos a [Eq. pU)) ]. Pro- 
vided that the slider 1 remains at rest, the force fl-^2 
decreases by = 2mg cosa(/is,2 — ^J■d) [by solv- 

ing the equation ([7])]. The slider 1 remains stable if 
^fi^2 + sina < ^s,i fng cos a [Eq. (|8])], i.e. if 
tana < /^d + (/^s,i — Ms,2)/2. We obtain a similar result 
by considering a contraction of the system. Thus, con- 
sidering the width of the distribution p(/^s) [Eq. pT|) ]. 
we estimate that the second slider generally remains sta- 
ble provided that tana < tana+ = 11^ — <7^/\/2, which 
defines a+ (Fig.El). 

The simplistic reasoning proposed above leads to the 
prediction of a well-defined plateau between a~ and a^ . 
However, a careful analysis of the numerical data reveals 
that < vc > slightly increases with the angle a even for 
large amplitudes Ag or ug (Fig.[21). To obtain the plateau 
values [Eq. (fT3|) ] or [Eq. p4)) ]. we implicitly assumed that 
any variation of 9 was enough for one of the conditions 
(IS]) and pUj) to be fulfilled, which is not correct when a 
dilation is followed by a contraction (or conversely). In 
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FIG. 4. Velocity < vq > vs. typical amplitude ag — In 

accordance with equation (|17|) . the average velocity < vq > 
tends to the asymptote /oo"s/(\/7rrth) only in the limit og 2> 
og (Here locjg is about 0.4 whereas the asymptote is reached 
to within 1% for l^ag about 100, Iq = 10^, /Xd=0.5, 7Zg=0.6, 
(7^ = 0.01 and tana = 0.25). 



this case, taking into account that the stable and unsta- 
ble sliders are swapped in Eqs. ([S]) and ([TU]). we obtain 



AXG 



The sliding dis- 



Zo|A6'| — 2(/igC0sa — sina) 

tance is reduced by the fact that a minimum dilation is 
necessary to invert the direction of the frictional force. 
As a consequence, in the case of periodic cycling, the 
velocity is given by: 



Tth < WG > — l^Ag — {fig cos a — sin a). 



(15) 



The relation ([T5)) holds valid provided that Ag > Ag 
defined by: 



Al = — (/ig cos a — sina). 
^0 



(16) 



If the amplitude is smaller than Ag, the dilations are 
not sufficient to rearrange the system and < vq >— 0. 
In the case of Gaussian variations [Eq. (IT2]) ]. a dilation 
(resp. contraction) is followed by a contraction (resp. 
dilation) in 2/3 of the timesteps. As a consequence, in a 
first approximation, the average velocity is given by 



Tth <VG > ' 



lo 



<7g — -{^J■s COS a — sina). 



(17) 



where we assumed that the dilations are enough to rear- 
range the system in each of the timesteps. The condition 
is reasonnably fulfilled when ag is much larger than: 



* 20F . 
ag = —j-[fi^cosa-sma). 

3*0 



(18) 



Thus, for ag ^ ag, < VG > increases linearly with ag 
[Eq. p7|) ]. When ag is deacreased, even below ag, the 
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dilations are always likely to rearrange the system and 
< vq > continuously decreases and vanishes for ag = 
(Fig.H. 



plateau velocities in the large amplitude limit. In ad- 
dition, we consider the dependence of the velocity on the 
amplitude of the temperature variations. 



B. Large systems: more than 2 sliders 



2. Analytic estimates 



Consider now a system consisting of N sliders con- 
nected by N-1 springs. It is important to notice that, for 
TV > 2, the system differs qualitatively from the system 
made of 2 sliders. Indeed, for N > 2, the dilation can 
induce the motion of the sliders at both ends without, 
necessarily inducing, a displacement of the slidcr(s) at 
center. Thus, the successive dilations do not necessarily 
induce a displacement of the center of mass in average. 



1. Numerical results 

In order to account for the creep of the system along 
the slope we consider the position of the center of mass 
xg = jj ^n=i ^ri as a function of time at the timesteps 
tq and report the average velocity < vq >, for large am- 
plitudes Ag or ag, as a function of the incline angle a 
(Fig. [S]). We observe that, for a small number N of slid- 
ers, < VG > exhibits a series of plateaus: < vg > takes an 
almost constant value in a finite range of the incline angle 
a. The number of plateaus increases when the number TV 
of sliders is increased. As expected, < vg > drastically 
increases when a approaches the critical angle ac- 
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FIG. 5. Velocity < vq > vs. incline angle a — The 

velocity < vg > exhibits a series of plateaus corresponding to 
constant velocities above critical angles onj. The number 
of plateaus increases with the number of sliders (dotted lines 
and open circles: N = 5; full lines, and full circles: A'^ = 6). 
The velocity < vq > diverges when a approaches the critical 
angle Oc (^o = 10^, fid=0.5, JX^—O.G, = 0.01 and ag =0.01). 

In the next section, we estimate theoretically the set 
of critical angles and the values of the corresponding 



The behavior of the system is understood by consid- 
ering the motion of the internal sliders. Let us assume 
that, during a dilation of large amplitude, a given slider i 
does not move whereas the sliders below move downwards 
and the sliders above move upwards. In the same way, 
let us assume that during a large contraction, a given 
slider j does not move whereas the sliders above move 
downwards and the sliders below move upwards. If the 
incline angle a is large enough, the sliders i and j differ 
(with i < j). In this case, the internal displacements re- 
sult in a reptation of the entire system along the incline 
as already described in the framework of the continuous 
models [i3,|Tl,[i3- 

Let us first determine the critical angle ; above 
which the slider i, previously at rest, starts moving down- 
wards, the slider i — 1 remaining at rest instead. Thus, 
let us first consider that, during a dilation, the slider i 
remains at rest. We assume, in a first approach, that 
the dilation rate is large enough for the sliders in mo- 
tion to slide continuously such that they are submitted 
to the dynamical frictional force ±^^"15 cos a. The con- 
dition is fulfilled provided that 9 ^ A/i/^oTdyn where we 
define A/i = fig ~ fJ-d- In addition, the dilation must 
be large enough for all the downstream sliders to move 
downwards. We shall later discuss this assumption. In 
this case, neglecting the inertia, one can indeed write the 
forces exerted by the slider i on the sliders i — 1 and i + 

fi^i-i + (* ^ 1) (Tng sin a + /id 'Tig cos a) = 
fi-yi+i + (TV — i){mg sina — jidfngcosa) = 

Considering the condition for the stability of the slider i, 

\fi-i^i + fi+i^i + m.g sin a| < fidmgcosa, 

we obtain that the slider i starts moving downwards 
above the critical angle ; given by: 



tan a 



1 



1 



TV 



(19) 



Note that j is a decreasing function of i such that the 
slider i — I remains stable for ; < a < 

Let us now consider a contraction of the same system. 
Assuming that for the chosen value of a, the slider i is 
at rest during the dilation, considering that the slider j 
remains at rest during the contraction, wc write: 

fj^j-i + [j — 1) (rngsina — ^d'mgcosa) = 
fj^j+i + (TV — j){mg sin a + ^^mg cos a) = 

Considering the stability of the slider j, replacing a by 
the critical value a^i' obtain that for a > j the 
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slider j = N + 1 ~ i starts moving downwards such that 
the sUder j + 1 is then at rest. 

In summary, for [ < a < q;Jj;_]^, the sHdcr i — 
1 remains stable during the dilation and the slider j + 

1 remains stable during the contraction of the system. 
Accordingly, we can write the displacements of any slider 
n, for a dilation Axn = [n - (i - 1)] loAe+ {A6+ > 0), 
for a contraction Axn = [n - {j + 1)] loA9^ {A9^ < 0). 
For a cycle A(?+ = —A9^ = 2Ag, one obtains the total 
displacement Ain = 2{N +3—2i)loAg for a total duration 

2 Tth. Note that Axn does not depend on n and that the 
associated velocity of the center of mass is simply 



rthwjf.i = {N + 3-2i)loAg. 



(20) 



We thus expect, for cycles of amplitude Ag, the plateau 
velocity ujj ; for ajj j < a < and 1 < i < {N + 3)/2 

(note that, for i = 1, ^ corresponds to the crit- 
ical angle of avalanche ac). For Gaussian tempera- 
ture variations Ag must be replaced by ag/^/n so that 
Tth^Nj = (A^ -|- 3 — 2i) l^agj^ in this case. We observe 
in the figure [5] that the equations and (PH)) give good 
estimates of the transitions and plateau velocities. 

When a contraction follows a dilation (or conversely) , a 
compression (resp. dilation) wave propagates from both 
ends inwards. Equation (^0)) is correct provided that the 
amplitude Ag is larger than a critical amplitude A*^ , such 
that the temperature variations induce the motion of all 
the sliders in the chain, which can be assessed in the 



following way: For ajj < a < a 
during a dilation, the slider i 
below downwards so that: 



and Ag > A*, 
1 pushes all the sliders 



h-i^t + [N-{i- 1)] (sin a - Hdcoaa) = (21) 

Note that during the dilation the sliders i to j are moving 
downwards so that the frictional force is already oriented 
upwards. During the contraction that follows, the sliders 
1 to j move downwards provided that the contraction 
wave propagating from the upper end reaches the sliders 
i — Let us now denote Ad* the corresponding dilation. 
For Ad*, the slider which is still at rest, pulls all 

the sliders above downwards so that: 

fi-i^i~2 + (i- 2)(sina - ^dcosa) = 0. 

The sliders i — 1 and i being still at rest, we get from 
equation (PT|) . 

ft^i-i = [iV - (i - l)](sina - ^dcosa) - IqAO* 

The temperature variation induces the motion of the en- 
tire chain provided that 

fi^i-i + f 1-2^1-1 + sin a - cosa > 

which leads to —IqAO* > A^(sinQ! — /id cos a) — A/icosa 
(remember here that Ad* < 0). Within a cycle, the 
amplitude of a contraction being A9 — —2Ag and the 



argument developped above holding true for a dilation 
following a contraction, we obtain the minimal amplitude 
leading to the reptation of the chain: 



Ag ~ {jid COS a — sin a) 
2h 



(22) 



where we assumed A^/N <C 1. Note that A'^ is inde- 
pendent of the position of the steady slider {i or j) and 
proportional to the number N of sliders in the chain. We 
can thus write, for finite amplitudes, 

~ \ N 1 
Tth^N,! = (^+3 — 2i) Zo ^ (/id cosa— sin a) . (23) 



The equation (j23p is in excellent agreement with the nu- 
merical results even when the static and dynamical fric- 
tional coefficients are not equal (A/i ^ 0, Fig.[S]). We thus 
also deduce from these observations that the average ve- 
locity is not sensitive to the width of the distribution 
of the static frictional coefficient. 




FIG. 6. Velocity < vg > vs. amplitude Aq or ag — 

For temperature cycles (open circles), < Vg > increases lin- 
early with the amplitude Ag above a critical amplitude A*g in 
agreement with the equation ^2'i\ [continuous line]. For ran- 
dom temperature variations (full circles), < Vg > increases 
significantly above a critical amplitude erg of the tempera- 
ture variations and reaches a linear asymptote for ag ^ erg in 
agreement with the equation (|24|) [continuous line. We remind 
that the asymptote is reached only when <jg ^ erg]. Random 
temperature variations thus lead to a smoother transition but 
do not change the qualitative behavior of the system (A^ = 30, 
lo = 10^, /id=0.5, /is=0.6, cr^ = 0.01, and tana = 0.25). 



In addition, the same qualitative behavior is expected, 
when the system is submitted to random temperature 
variations. In the figure [SI we observe that equation ([25)1 
agrees with the numerical results, provided that the am- 
plitude Ag is replaced by ag j ^/n and Ag by <7g = ^y/nAg 
so that: 



rth«Ni = (^+3-2z)[o 



ag N 

—= ^ (/id cosa— sin a) (24) 

3^0 
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from equation ((20)) . provided that fc = NYS/L, 



\ 



100 



150 



200 



FIG. 7. Velocity < va > vs. number of sliders TV — 

For temperature cycles (open circles) and random tempera- 
ture variations (full circles), < Vg > exhibits a non-monotonic 
behavior as a function of A''. For small N , < Vg > increases 
because of the increase of the size of the system; for large A'', 

< Vg > decreases because of the increase of critical amplitude 
Ag or ag. For temperature cycles, the numerical results are in 
excellent agreement with the equation (|23|l (continuous line). 
Note however the dispersion of the numerical points for small 
N which is due to the fact that, at a given amplitude As, 

< Vg > corresponds to discrete plateau values which are not 
accounted for by the equation (|23|l when considering contin- 
uous values of a {lo = 10"^, /id =0.5, 7ij,=0.6, cr^ = 0.01, Ag — 
0.005 or ae = 0.005^^ and tana = 0.25). 



for random temperature variations. 



IV. DISCUSSION 

In the continuous description proposed by Moseley and 
Bouasse [13 j the elasticity of the material was ne- 
glected which led to the conclusion that the creep veloc- 
ity was proportional to the amplitude of the temperature 
variations. The introduction of the elastic effects leads 
to the conclusion that the temperature variations induce 
the creep of the system only if their amplitude is large 
enough, as already proposed by Croll. However, note 
that the critical value given by the equation differs 
from the result proposed in The difference comes 

from the fact that Croll considered that the system was 
free of stress at the beginning of each phase of a cycle 
(dilation or contraction). In our approach, each dilation 
(resp. contraction) follows a contraction (resp. dilation) 
and the frictional force is initially non zero, mobilized 
in the opposite direction. Considering that our chain of 
sliders model a solid of mass M, length L, cross section 
S made of a material having a Young modulus Y and 
a density p creeping along an incline, one can estimate 



TthVG 



L- 



tana 



~ (,Md COS a - sm a ) 



(25) 



We thus obtain that the creep velocity is independent of 
the number TV of contacts. The minimum amplitude of 
the temperature changes that produce the creep of the 
system along the incline can be written: 



AT* = — — {ula cos a — sin a) 
TYk ' 



(26) 



Thus, for a given i, the system is thus more likely to 
creep along the incline for a larger Young modulus Y , a 
larger thermal expansion coefScient k and a larger angle 
a. Note again that the stability of the system does not 
depend on the number TV of contacts. Typically, with 
i = 1 cm, y - 100 GPa and p ~ 10"^ kg.m'^, we es- 
timate that relative dilations of about 10~® are enough 
to make the system creep. Such dilations, which corre- 
spond to temperature changes of about AT ^ 1 niK, are 
difficult to avoid and, in general, any frictional contact 
between two macroscopic solids cannot be considered as 
perfectly static. 




500 



FIG. 8. Position < xq > vs. time t — For Ag ~ A%, the 

system exhibits an irregular dynamics even when temperature 
cycles are imposed. Note that fast creep can be followed by 
long quiescent periods as observed for Aq = l.OlAg (black) 
and Ag = 0.99 (grey) for A% ~ 3.6 lO"'' (TV = 30, h = 10^ 
/id=0.5, /is=0.6, cr^j = and tana = 0.25). 

The second practical situation mcntionncd in the intro- 
duction is the creep of g ranular material induced by tem- 
perature changes |7Hl3j. In this case, the length would 
account for the typical grain size, k for the rigidity asso- 
ciated with the grain-grain contact and TV for the typical 
number of grains in one typical dimension L = TV^o of the 
system. It is then particularly interesting to consider the 
dependence of the typical creep velocity on the number TV 
[Eq. , Fig. [7] . For a given amplitude of the tempera- 
ture variations, < Vg > exhibits a non-monotonic behav- 
ior as a function of TV, increasing linearly with TV for small 
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N and decreasing for large N because of the increase of 
the critical amplitude of the temperature variations likely 
to induce the creep [Eq. ()22|) ]. For temperature cycles, 
from equation (j22[) . < Wg > is expected to vanish above 
for TV > 2lQAg/ (fid cos a — sin a) whereas random tem- 
perature variations are expected to be always likely to 
produce creep. 

Let us now focus on the dynamics of the system sub- 
jected to temperature cycles of amplitude close to Ag 
(Fig. [8|). We observe that, for Ag close to Ag (above 
and below Ag) the system creeps in an irregular man- 
ner. The position of the center of mass, ic, exhibits a 
series of rapid variations (jumps), separated by periods 
of time during which the system is apparently at rest. 
This latter conclusion holds true even if the value of the 
static frictional coefficient is well-defined, i.e. even for 
CT^ = 0. Thus, for Ag approaching Ag, one observes a 
transition from a continuous creep regime, during which 
XQ decreases at each cycle, to the irregular creep regime. 
The transition is reminiscent of the transition between 
a continuous-flow and an intermittent-flow regime, ob- 
served when the amplitude of the temperature cycles im- 
posed to a granular column is decreased [l0l - [l3| . The ex- 
perimental and theoretical systems are very different but 
there are some common aspects like the frictional contact 
between the particles that are in a limited number N. We 
can however attempt to estimate the critical amplitude 
Ag expected from the model for a column (diameter 1 
cm) of glass beads (typically 500 fim in diameter). The 
elasticity of the system is due to the Hertz contact be- 
tween the grains and we can estimate that the stiffness 
k depends on the pressure. Denoting S the penetration 
distance and R the radius of the grains, we can write 
k ^ Y{RSy^'^. For an infinite vertical column, because 
of Janssen effect, one can estimate the local pressure 
P ^ pgD, where D stands for the diameter of the column 
and p for the density of glass. Writing that the force ap- 
plied to the grains kS ~ PR'^, we get k ^ R[pgDY'^]^/^ . 
From equation with Iq 2R and N D/2R, we 

estimate A^ ~ {pgD /Y)"^/^ . Note first that the result 
does not depend on the size of the grains and, thus, not 
on the number of grains in the diameter of the column. 
Second, with p^2 10^, D - 1 cm and F 20 GPa, we 
obtain that - 5 10"^ and, thus, with k ~ 3 10"^ K^^ 



for glass, that a transition between the continuous- and 
the irregular-flow regimes is expected for amplitudes of 
the temperature variations of the order of 1 K. The lat- 
ter value is interestinglyclose to the experimental value 
which is of about 3 K [10[ . Even if the good agreement 
between the theoretical estimate and the experimental 
values might be accidental, we think that the compari- 
son between our model and the granular column is worth 
to be mentionncd. 

V. CONCLUSION 

We reported on the detailed behavior of a frictional 
system subjected to thermal dilations. We observed that 
for a small number of sliders in the chain and large am- 
plitude of the dilations, the 'reptation' velocity of the 
center of mass exhibits a series of plateaus as a function 
of the incline angle. In the limit of an infinite number 
of sliders and large amplitudes, we recover former results 
obtained for the creep of a solid on an incline. Because 
of the elasticity of the material, the creep velocity is ex- 
pected to vanish for a finite value when the amplitude of 
the cycles is decreased. We obtain an expression of the 
critical velocity which slighly differs from former results. 

Finally, for a finite number of sliders, we observe nu- 
merically that the system experiences an irregular trajec- 
tory, the center of mass sliding rapidly between quiescent, 
rather long, periods of time even if the system is sub- 
jected to periodic cycling. The transition between the 
continuous and the irregular creep depends on the size 
of the system. The irregular creep is reminiscent of re- 
cent observations of the irregular compaction of granular 
matter under the action of periodic temperature changes. 
The systems are very different but we do believe that the 
study of this specific regime will provide us with interest- 
ing clues for the understanding of the peculiar behavior 
of granular matter. The latter study will be the subject 
of a further publication. 
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